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Abstract 



The work described in this paper is the first step toward a relativistic 
three-quark bound-state calculation using a Hamiltonian consistent with the 
Wigner-Bargmann theorem and macroscopic locality. We give an explicit 
demonstration that we can solve the two-body problem in momentum space 
with spin-dependent interactions. The form of the potential is a combina- 
tion of linear+Coulomb+spin-spin-|-spin-orbit-|-tensor, which includes con- 
finement and is of the general form consistent with rotation, space-reflection 
and time-reversal invariance. Comparison is made with previous calculations 
using an alternate technique and with the experimental meson mass spectrum. 
The results obtained suggest that the model is realistic enough to provide a 
two-body basis for the three quark baryon problem in which the Poincare 
group representation is unitary and cluster separability is respected. 

PACS number(s): 12.40. Qq, 12.40.Aa 
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I. INTRODUCTION 



The work described in this paper is the first step toward a relativistic three-fermion calcu- 
lation. For light quark systems, wherein mesons are described by quark- ant iquark dynamics 
and baryons by three quark dynamics, one calculable approach is the relativistic potential 
model. Covariant approaches based on the Bethe-Salpeter equation |J appear to require 
some form of reduction to a soluble equation, e.g., Salpeter's instantaneous approximation 
0, ladder approximations in Euclidean space pj or a relativistic equation in which one of the 
quarks is restricted to its mass shell [Q. Each of these reductions can be argued to have its 
own strengths and weaknesses. There is currently no truly satisfactory solution of the Bethe- 
Salpeter equation for bound states of two fermions, let alone three. Therefore, we resort to 
a potential model with a model Hamiltonian consistent with the principle of relativity and 
macroscopic locality || to solve the three light-quark bound-state problem. The three-body 
mass operator (Hamiltonian) can be defined using the Bakamjian- Thomas construction ||. 
Unfortunately such a relativistic Hamiltonian is already very difficult to solve 0,0. One 
feasible approach is to solve the Faddeev equations in momentum space using appropriate 
two-body wave functions as a basis [0,|8|]. This necessarily requires an accurate calculation 
of such a two-body wave function basis. We have chosen the collocation method to solve 
the two-body equation rather than variational techniques for this reason. Consequently, 
it is necessary to first solve the quark-antiquark problem accurately in momentum space 
(because of the difficulty arising from the square root operator in coordinate space using the 
collocation method) with general spin-dependent interactions. This has not previously been 
done. Previous work in momentum space has been limited to spin-independent calculation 
|7]U. In the present work we attempt to ensure that the two-body basis is optimal, we fit the 
meson masses (i.e., quark-antiquark systems) with a parametrized interaction in momentum 
space including a full range of spin-dependent interactions. 

In section II we describe the relativistic two-body Hamiltonian in terms of a linear con- 
finement term , a Coulomb term, and various spin-dependent pieces. We emphasize that this 
is not a simple nonrelativistic reduction of an effective one-gluon-exchange potential since 
the coefficients of the various terms in the potential will be constrained by phenomenological 
considerations alone. The various terms making up the potential are used because they are 
invariant under rotations, space-reflection and time-reversal. To treat singularities in the 
spin-dependent interactions a form factor is introduced for the quark-gluon vertex. Since 
the nonperturbative quark-gluon vertex is in general momentum-dependent, this is a natu- 
ral thing to include. Numerical solutions are obtained for the model of reference |I(J using 
partial wave expansions of the various interactions given in section II. Results for various 
parameter choices including those of reference |Ttl are presented in section III. These re- 



sults show that the light meson spectrum can be well described by such a model. Planned 
future directions involving other properties of the meson system, such as electromagnetic 
form factors, and the extension to relativistic three-fermion systems are given in section IV. 



II. FORMALISM 

In our relativistic potential model, the Hamiltonian H is the sum of a relativistic kinetic 
energy operator T and a potential operator V. The kinetic energy operator has the form 



2 



T = J Ml + k 2 + 



Mf + k 2 , 



(1) 



where Mi, Mi are the masses of the two quarks and k is the relative momentum. In what 
follows, we assume M = M\ = M2. The potential is the sum of a linear confining potential 
and a short range spin-dependent interaction: 



V = V L + V S . 



(2) 



Here we use a linear confining potential Vl = \for where \fa is the string tension and r 
is the relative coordinate. As a guide to choosing an effective spin-dependent short range 
interaction we begin by noting that the nonrelativistic reduction of the one-gluon-exchange 



potential in momentum space has the form [|TT 
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where a s is the strong- interaction fine-structure constant, f c is the color factor (which is 
—4/3 for quark-antiquark and —2/3 for quark-quark), cr^ cr 2 are the Pauli matrices and 
q = k' — k is the momentum transfer. It is readily seen that the four terms in the curly 
brackets represent the Coulomb potential, the spin-spin, spin-orbit and tensor interactions 
respectively. These are easily recognized in coordinate space representation: 
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where S = |(cn + er 2 ) is the total spin operator and L is the orbital angular momentum 
operator. The spin-spin interaction includes a delta function of r, (i.e., it is a contact inter- 
action), and so we introduce a Gaussian form factor exp(— |A 2 q 2 ) at the quark-gluon vertex 
as in Ref. pTJfl . The variable A can be interpreted as the size of the quark. In Ref. [T2| 
a form factor l/(q 2 + (3 2 ) in which is the effective quark size is used to eliminate the 
singularity. Since the one-gluon-exchange potential (Eq. ||]) is derived via a nonrelativistic 
reduction (and so can not represent a reasonable interaction in a relativistic calculation) 
and in any case is really an effective interaction, we are completely free to vary the relative 
strengths of these interactions. In particular it is well known that the phenomenological 
strength of the spin-orbit interaction is much weaker than that predicted by the nonrela- 
tivistic reduction of the one-gluon-exchange potential . Hence we use multiplying factors 
Cl, Cq, Css, Cls, Ct so that the strength of every term can be varied to fit the meson data. 
The general form of the potential that we used in our calculation then has the form 
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(5) 



Note that the only potentials that respect rotation invariant, space-reflection and time- 
reversal which are not included in our potential are ((T\ ■ L)(cr 2 • L) and (<Xi — cr 2 ) • L. The 
latter term becomes significant only if I ^ and one of the quarks is much heavier than 
other. 

To do the calculations in momentum space we have to use a partial-wave expansion 
for the potentials in Eq. (|5]). The difficulty with a momentum-space formulation for the 
linear potential was solved in Refs. |J14] , |l5|| . The following are the results for the partial- wave 
expansion of the potentials in momentum space. 
For potentials depending on |q| only: 
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For the spin-spin interaction: 
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For the spin-orbit interaction: 
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and for the tensor interaction: 




{ J 2 n}EEHf +3+1 (2^i) 
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is a 6-j symbol [0. Note that equa- 



In the above equations CJ} is a binomial coefficient, Pi(fi) is a Legendre polynomial, 

( jl 32 3: '' ) is a 3-j symbol and { ' 32 33 ' 

tions (8b~) and (95) are new results. With these two equations, we can do calculations with 
spin-orbit and tensor interactions for almost any functional form of V(|q|) in equations 
and (H). 

To solve the relativistic wave equation 
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we first introduce an auxiliary function ipi S j(k) given by 



(10) 



(11) 



The power / — 1 follows from considering the behavior of ^i S j(k) as k —* 0. The auxiliary 
wave function is expanded in a complete set of basis functions. In Refs. Jl2| and [l3j the 



wave functions were expanded in a harmonic oscillator basis. In Ref. [pj]] the wave functions 
were expanded in confinement eigenstates which are solutions of a relativistic Schrodinger 
equation with a confinement potential only. Here we choose cubic Hermite splines 



as 



the basis functions. The cubic Hermite splines are piecewise polynomials of degree three 
with continuous first derivatives [18]. Also we impose the boundary conditions by requiring 
the auxiliary wave function be zero at k = and at a cutoff k = k max . By expanding the 
wave function in N splines functions and requiring that the expansion satisfies Eq. (|TUp at N 
distinct values (collocation points) of k Eq. (pi]) can be converted to a matrix equation. The 
two-point Gaussian quadrature points on each interval are chosen as the collocation points. 
The eigenvalues and eigenvectors of the matrix equation are the masses of the mesons and 
the coefficients of the meson wave function bases respectively. 



III. RESULTS 

The light meson mass spectra can be obtained by solving for the eigenvalues of the 
relativistic Hamiltonian H = T + V — M qq - where V is given in Eq. @ and M qq - is a constant 
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energy to be determined from data. If M q q = 0, the model gives the meson mass spectra; 
if Mqq ^ 0, then the model only gives the splittings. In our opinion, this overall constant 
energy is somewhat artificial. In our first calculation, the parameters used are M = 360 
MeV, f c a s = —0.5 and ^ = 0.197 (GeV) 2 . These parameters are chosen so that the results 
can be compared with the results of Ref. [[TIJ which used a variational method to solve the 
wave equation in coordinate space. Unfortunately, they used a different method to treat 
the spin-orbit interaction so we can only directly compare our results with theirs without 
spin-orbit interaction, i.e. for n,p and b mesons. We initially set Cl = Cc = Css = 1, 
Cls = Ct = and eliminate the vertex form factor for the Coulomb potential. These 
results are listed in Table | with the same constant energy (M qq - = 750 MeV) subtracted 
out. 



Our results are consistently higher than their results by 25 — 60 MeV |I9| . We have been 
unable to establish the origin of the discrepancy but note that our Fortran codes were cross 
checked with a totally different method |2D| and we agree within numerical error (i.e., < 0.5 
MeV). 

The results of fitting light meson spectra by varying two parameters A and Cls (Cl = 
Cc = Css = Ct = 1, M = 360 MeV, no form factor for the Coulomb potential) are shown 
in Table |TJ. The differences between this second calculation and that of Ref. |JTU] are due to 
the following factors: 



a) We use A = 0.6426 GeV _i ; they used A = 0.13 fm = 0.658805 GeV . 

b) The spin-orbit interaction is treated differently and our strength of the spin-orbit 
interaction is Cls = 0.3236. 

c) We consider the coupling between different channels via the tensor interaction. 

d) A different constant energy (M qq - = 738.24 MeV) is used here versus 750 MeV. 

The average deviation from the fifteen experimental values is 69 MeV for our results and 92 



MeV for Ref. Pf 



As we mentioned in the previous section, we are free to adjust all seven parameters (M, 
A, Cl, Cc ,Css, Cls, Ct) to fit the light meson spectra. We also eliminate the need for 
the constant energy term M qq used above and use a form factor not only for spin-spin, spin- 
orbit and tensor interaction but also for the Coulomb interaction. The parameters we use are 
M = 258 MeV, A = 0.645 GeV" 1 , C L = 0.6704, C c = 3.3824, C ss = 0.155, C LS = 0.0448 
and Ct = 0.1868. The results of this fitting are listed in the last column of Table [TI|. 
The average deviation is 28 MeV for eleven states which can be compared to the average 
deviation of 55 MeV in a similar model calculated in Ref. ]T2|] . The average deviation from 



all fifteen experimental values is 51 MeV which is significantly better than that obtained by 
only varying two parameters (69 MeV). 

The above parameter values reflect the following: The linear confinement potential is a 
little weaker than for static quarks (Cl — 1). The Cc value is considerably greater than 
one, which means that light quark systems have a running coupling constant (a s Cc) much 
larger than heavy quark systems as might be expected from QCD (asymptotic freedom). 
Cls is very weak as observed by others [ IB] . Css and Ct are also reduced relative to values 
obtained for heavy quarks. 
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IV. FUTURE WORK 



Potential models based on the constituent quark model have enjoyed considerable success. 
Even the nonrelativistic potential model works much better than one might naively expect. 
However, the nonrelativistic constituent quark model fails to explain the meson and baryon 
form factors at high energy because the calculated form factors at high energy fall off too 
fast. In our calculations, the tails of the momentum space wave functions have a power-law 
fall off, rather than the exponential decay of harmonic oscillator potential models. It is 
hoped that using a wave function like ours with an appropriate relativistic "boost" might be 
able to provide the correct form factor behavior at high energy. This question is currently 
being pursued. 

In 1939, Wigner showed that the description of a physical system in a relativistic 
quantum theory must correspond to a unitary representation of the Poincare group. In 
1961, Foldy recognized the importance of macroscopic locality (cluster separability) as 
an additional constraint on relativistic potential models. Many "relativistic" three-body 
potential models |23|] do not satisfy either of these two requirements. Relativistic three- 
body mass operators which satisfy both of these two requirements have been discussed in 
a recent review ||. A relativistic three-body bound state formulation which incorporates 
Wigner's theorem and macroscopic locality can be solved by using the relativistic Faddeev 
equation for the eigenvalues and eigenvectors of the three-body mass operator in momentum 
space. The calculations are very complicated since the potential operators are functions of 
the position operators and are inside the square roots. These difficulties can be overcome 
by using relativistic two-body wave functions as the basis J7||§. Therefore, solving the two- 
body equation accurately in momentum space was necessary as a first step. The two-body 
calculations discussed in this paper are important primarily because of their use as input 
for three-quark calculations. Our ultimate purpose is to solve the relativistic bound state 
problem for three light constituent quarks. 
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TABLE I. The light meson mass spectrum without spin-orbit and tensor interactions. Our 
results (the fourth column) are consistently higher than those of Ref. JE0] by 25-60 MeV. 
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TABLE II. The light meson mass spectrum by varying two parameters A and Cls (t ne fifth 
column) and those by varying all 7 parameters (the last column). The average deviation from the 
fifteen experiment values is 69 MeV for the second calculation (the fifth column) and 51 MeV for 
the third calculation (the last column). 
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